Fiber-optic realization of anisotropic depolarizing quantum channels 
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We employed an electrically-driven polarization controller to implement anisotropic depolarizing 
quantum channels for the polarization state of single photons. The channels were characterized by 
means of ancilla-assisted quantum process tomography using polarization-entangled photons gener- 
ated in the process of spontaneous parametric down-conversion. The demonstrated depolarization 
method offers good repeatability, low cost, and compatibility with fiber-optic setups, ft does not 
perturb the modal structure of single photons, and therefore can be used to verify experimentally 
protocols for managing decoherence effects based on multiphoton interference. 
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I. INTRODUCTION 

The practical success of quantum technologies relies on 
methods to protect and to enhance quantum effects in the 
presence of decoherence. For example, imperfect entan- 
glement between distant parties can be improved using 
rotocols for entanglement purification and distillation 
l| . Transmission and storage of quantum superpositions 
can be made more robust with the help of quantum error 
correction techniques 0] . Demonstrating experimentally 
the feasibility of these methods is a crucial step in the 
development of quantum technologies. 

In this paper we present an experimental realization of 
a depolarizing fiber optical channel that can be used to 
simulate a range of non-trivial decoherence mechanisms. 
We anticipate that the constructed setup can provide a 
testbed to verify a number of theoretical proposals for 
managing decoherence effects. As optical methods for 
generating and manipulating entangled states are rela- 
tively well developed, photonic systems offer a promising 
route to implement more complex protocols for protect- 
ing quantum information in the presence of noise. Opti- 
cal approach has been already used to verify experimen- 
tally decoherence-free subspaces 0, Q , use them in ele- 
mentary quantum computing Q and quantum cryptog- 
raphy [6] protocols, and to carry out single-qubit purifi- 
cation |7i and quantum error correction [Sj]. The device 
presented in this paper can be programmed to simulate 
specific anisotropic depolarizing channels that appear in 
practical situations, for which more efficient specialized 
protection methods may exist [l^ . 

One practical method to generate decoherence in opti- 
cal qubits is to introduce correlations between the mea- 
sured degree of freedom — typically polarization — and an- 
other one, such as the frequency [3] or the wave vec- 
tor [ll|, that is not sensed by the measuring apparatus, 
and from the formal point of view can be considered as 
the environment. This approach however modifies the 
modal structure of the fields. A polarization qubit pre- 
pared initially as a superposition of one photon in two 
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modes that differ only by their polarization and are iden- 
tical otherwise, becomes a complex object that cannot 
be assigned a single spatio-temporal profile. This has 
profound implications when such qubits are manipulated 
using presently the most accessible experimentally tool 
of multiphoton interference realized in linear-optics net- 
works. Such interference requires excellent matching of 
the spatio-temporal modes, and its lack provides distin- 
guishing information that deteriorates interference visi- 
bility Iisj . For demonstrating more complex manipu- 
lations of noise-affected states based on multiphoton in- 
terference, quantum channels in the polarization degree 
of freedom need to be implemented without changing the 
spatio-temporal properties of the fields, and the approach 
presented in this paper satisfies this requirement. In our 
case, the effective quantum channel is obtained through a 
temporal average over suitably chosen polarization trans- 
formations. This approach simulates decoherence rather 
than induces it through an interaction with another mi- 
croscopic ancilla system, with a completely equivalent 
outcome. 

Our realization of a depolarizing channel is a fiber op- 
tic counterpart of the free-space scheme used in 01 based 
on Pockels cells. Transformations of the state of polariza- 
tion (SOP) are introduced with the help of an electrically 
driven polarization controller composed of magnetic ele- 
ments that introduce birefringence by squeezing mechan- 
ically the fiber • The controller is an inexpensive low- 
voltage device, compatible with all-fiber setups and char- 
acterized by negligible insertion losses. We verify the ac- 
tion of the channel through ancilla-assisted quantum pro- 
cess tomography using maximally entangled photon 
pairs generated in the process of spontaneous parametric 
down-conversion [iB]. Sending one photon through the 
channel and performing quantum state tomography on 
the entire two-photon state yields complete characteriza- 
tion of the channel thanks to the Jamiolkowski isomor- 
phism |171] . 

We note that scrambling the SOP of optical fields is 
an important technique in fiber optic communications, 
which can reduce the impact of deleterious polarization- 

tendent effects in practical transmission systems [isl - 
Compared to applications in high bit rate fiber 
communication systems, requirements set by proof-of- 
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principle demonstrations of quantum protocols are some- 
what different. Firstly, shortening time scale of polariza- 
tion scrambling is not of primary importance, as typ- 
ical experiments with entangled photons require long 
collection times. Therefore mechanical scramblers are 
presently well suited for this class of applications. Sec- 
ondly, the actual distribution of the SOP after scrambling 
is not important, as only ensemble averages matter for 
experiments in which consecutive outputs are manipu- 
lated and detected individually. 

The structure of this paper is as follows. In Sec. Ullwe 
summarize the theoretical basis of our experiment. The 
experimental setup is presented in Sec. IIIH and obtained 
results are discussed in SecHV] Finally, Sec. IVl concludes 
the paper. 



II. THEORY 

Let us start with a brief theoretical review of qubit 
channels. The most convenient description is based on 
the Bloch vector representation. A general channel £ 
is given by an affine map acting on the Bloch vector r 
according to plj : 




FIG. 1. The tetrahedron of singular values for unital quantum 
channels which results from the physical condition of complete 
positivity. Dots and squares represent channels implemented 
and characterized experimentally by means of ancilla-assisted 
quantum process tomography, discussed in detail in Sec. IIVI 
Roman labels (a-d) refer to channels depicted in Fig. (3] 



ri-^Ar + a, (1) 

where A is a real 3x3 matrix and a is a vector character- 
izing the output state for a completely mixed input. We 
will be interested here only in unital maps, which do not 
change the maximally mixed state, and therefore a = 0. 
Singular value decomposition enables us to represent A as 
a product A — ODO' of two proper rotations O, O' and 
a diagonal matrix D = Aia,g{Dx,Dy,Dz). This allows 
us to separate the trivial rotations of the Bloch vector 
corresponding to unitary transformations from the depo- 
larization process, described by the three singular values 
comprising D. Distinctly from the usual convention we 
allow here negative singular values, which stems from the 
requirement that the rotations O, O' are proper. 

Physical channels must satisfy the condition of com- 
plete positivity, which means that a trivial extension of 
the channel £ to any auxiliary system yields a positive 
map. For unital maps this condition implies the following 
constraints on D [21]: 

\D,±Dy\<\l±D,l (2) 

where the ± signs on the left- and the right-hand side 
of the above equation are taken independently. Geomet- 
rically, the set of allowed parameters D^, Dy, Dz forms 
a tetrahedron shown in Fig. [TJ whose vertices are given 
by the identity I and 180° rotations 'Rx,'R-y,'Rz about 
axes of the coordinate system in the Bloch representa- 
tion. These rotations correspond to unitary transforma- 
tions described by Pauli matrices ax-,^y,^z- The geo- 
metrical picture shows that the non-unitary depolarizing 
part of an arbitrary unital channel given by D can be 



represented as a convex combination of the identity and 
three rotations: 

D = Pol + Px^x + Py'R.y + Pz'R-z, (3) 

where the nonnegative parameters po,Px,Py,Pz form a 
probability distribution. The decomposition given in 
Eq. Q defines a practical way to implement an arbi- 
trary unital channel by a random applications of identity 
I and rotations Ra:,Rj,,R2 with appropriate probabil- 
ities, sandwiched between unitary transformations that 
realize O and O'. This approach will be the basis of our 
experiment. 

We will characterize the action of the channel using 
ancilla-assisted quantum process tomography, which is 
based on the Jamiolkowski isomorphism [l7| between 
quantum channels and bipartite states. Specifically, 
preparing two qubits in a maximally entangled bipar- 
tite state |<I>) and sending one of the subsystems through 
the channel yields a state gs = (£ (8) /)(|$) ($|), that car- 
ries complete information about the channel. Performing 
quantum state tomography on gg allows one to charac- 
terize fully the channel £. The density matrix of the 
maximally entangled input state can be written in the 
form [13: 

|1>)(<I>|-Mi®i-5]T,,<T,®aJ , (4) 

where T is a real 3x3 matrix with the singular 
values equal to ±1. For example, a state \^+) = 
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FIG. 2. Experimental setup. BC, Babinet compensator; 
BBO, a pair of beta-barium borate crystals; IF, interference 
filter; F, red filter; MFC, manual fiber polarization controller; 
EPC, electrically-driven fiber polarization controller; PBS, 
polarizing beam splitter; Dl, D2, single-photon detectors. 



lt^))/V2, close to the one used in our exper- 
iment, corresponds to T = diag(— 1, 1, — 1). For unital 
channels, the output state qs has a form analogous to 
Eq. with the matrix T replaced by AT. Its ele- 
ments can be calculated as expectation values (AT)ij = 
Tr[^£(iJi (Tj )]. Multiplying the reconstructed matrix by 
yields A. This method tolerates imperfections in 
the input state such as mixedness, as long as the matrix 
T is invertible. Because our primary interest here are 
depolarizing properties of quantum channels and unitary 
rotations can be easily compensated with polarization 
controllers, we will use experimental data to extract the 
singular values of A that are responsible for the non- 
unitary dynamics. 

It is worth noting that applying the depolarization 
method described above to a maximally entangled state 
allows one to generate an arbitrary state of the form given 
by the right hand side of Eq. Q . The matrix T can be 
brought to the diagonal form analogously to the Bloch 
representation of unital channels, with the singular val- 
ues satisfying the analog of Eq. ^ which stems from 
the positivity condition for density matrices. The output 
state can be fully controlled through voltages applied to 
a single device, unlike sources of mixed bipartite states 
which rely on more complex preparation procedures (23| . 

III. EXPERIMENTAL SETUP 

The experimental setup is shown in Fig. [2j The ultra- 
violet pump beam entering the setup from the left is pre- 
pared by doubling femtosecond pulses from a titanium- 
sapphire laser (40 fs FWHM pulse duration, 80 MHz rep- 
etition rate, 780 nm central wavelength) in a beta-barium 
borate (BBO) crystal to yield second harmonic with the 
average power of 14 mW, which is separated from the 
fundamental with two dichroic mirrors and a blue color 
filter. 

The source of entangled photon pairs is based on spon- 
taneous parametric down-conversion in a pair of type-I 
BBO crystals oriented at 90° with respect to each other 
[l6l] . The half-wave plate A/2 and the Babinet compen- 



sator BC preceding the crystals ensure that the pairs 
produced by both crystals are temporally indistinguish- 
able and have equal probability amplitudes. The gener- 
ated signal and idler photons are sent through red filters 
F and 5 nm bandwidth interference filters IF to filter 
out the scattered pump field and to restrict the spectral 
bandwidth. 

The idler photons are then coupled into a single-mode 
fiber using an aspheric lens. They pass first through a 
manual fiber polarization controller MFC with settings 
fixed at the beginning of the experiment to compensate 
for the overall SOP change in the fiber, and then through 
the electrically driven polarization controller EPC (OZ- 
Optics model EPC-400) which implements the actual de- 
polarization channel. The idler photons are subsequently 
coupled out of the fiber into a collimated beam with an 
aspheric lens and sent through a quarter-, a half- wave 
plate and a polarizer to select a specific polarization 
state, detected finally with a free-space single-photon 
counting module Dl (Perkin Elmer model SPCM-AQ- 
131). An identical set of two wave-plates and a polarizer 
is placed in the path of the signal photons before coupling 
them into a single-mode fiber connected directly to a 
module D2 (Perkin Elmer model SPCM-AQR-14-FC). To 
reconstruct the density matrix, coincidence count rates 
were measured in a basis of 16 different two-photon po- 
larization states that form a tomographically complete 
set. In the standard notation, the basis consisted of the 
following states: \^X), \XX), |0o), |0t), 

l^t), 1^^), 1^0), 1^^), 10^), it^), ito), 

1^0), |t>0), selected by appropriate settings of quarter- 
and half-wave plates in both arms. The tomographic re- 
construction of the density matrices was carried out using 
the maximum likelihood technique [13, HE] • 

The electrically driven fiber polarization controller 
consists of four magnetically driven squeezers driven by 
external voltages, which apply mechanical stress to the 
optical fiber thus inducing birefringence. Each squeezer 
works as a wave-plate with a fixed orientation and a tun- 
able phase delay. It is easy to see that placing a sequence 
of three such squeezers with relative 45° rotations be- 
tween their optical axes allows one to implement an ar- 
bitrary rotation of the Bloch sphere, with phase delays 
defining the three Euler angles of the rotation [l3| . The 
device used in the experiment included four squeezing 
elements, and we found that their relative orientation 
required more complex driving to achieve necessary po- 
larization rotations. Switching between several different 
settings of applied voltages and averaging over time yields 
a depolarizing channel with controllable characteristics. 
The time of averaging was determined by the brightness 
of our source of entangled photon pairs and was chosen 
to be 70 s for each setting of the quarter- and half-wave 
plates needed for density matrix reconstruction. That 
time was much larger compared to the duration of a sin- 
gle setting of the polarization controller, which was of the 
order of 100 ms (exact figures will be specified in Sec. lIVI) . 
By employing a polarization controller faster then the one 
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used in our experiment |2Q|, this depolarization method 
could be used with much shorter acquisition times and 
thus with brighter sources. The key advantage of this 
depolarization method is that it does not introduce cor- 
relations between the polarization state and spatial or 
spectral degrees of freedom of transmitted photons, as 
long as polarization mode dispersion can be neglected 
over the spectral range of the transmitted photons. 

We define the quantum channel whose characteristics 
are to be measured as the polarization transformation 
experienced by the idler photons between the aspheric 
lenses coupling into and out of the fiber. The action 
of EPC was calibrated by coupling into the channel a 
macroscopic laser beam prepared in three different states 
of polarization: horizontal, diagonal, or left-circular, and 
performing a standard SOP measurement on the out- 
put using an ellipsometer that detected intensities of 
horizontal- vertical, diagonal, and circular components on 
three parts split off the input beam. In the following, we 
will refer to this measurement as simple process tomogra- 
phy to differentiate it from the ancilla-assisted quantum 
process tomography, that is the core of our experiment. 
In the first step of the calibration procedure, the driving 
of the EPC was switched off completely and MPC was 
aligned manually such that the action of the channel was 
equal to the identity. Starting from this reference state, 
voltage settings required to achieve 180° rotations were 
determined. 



IV. RESULTS 

The input state for ancilla-assisted quantum process 
tomography was prepared by setting the action of the 
quantum channel to identity and measuring the density 
matrix of the generated state. For measurements de- 
scribed in this paper, coincidence counts were collected 
over 70 s for each pair of polarization settings. The re- 
sults of quantum state tomography performed on the in- 
put state are shown in Fig. ^a.). The linearized entropy 
of this state is equal to 0.1 ± 0.1, and the amount of en- 
tanglement can be quantified with the concurrence value 
of 0.92 ± 0.07. The input state is well approximated by a 
pure state of the form |$+) = -hexp(i(^) \tt))/V^ 

with the phase equal to <j) = —0.12. The fidelity between 
the measured state and its pure-state approximation is 
equal to 0.98 ± 0.09. 

We concentrated our efforts on the implementations of 
anisotropic channels. In the first series of measurements 
we studied channels that lie on the edge of the tetrahe- 
dron, and can be realized as a mixture of the identity I 
and a single 180° rotation R, which is obtained by ap- 
plying an appropriate voltage to one of the squeezers in 
the EPC. Thus, the realized quantum channel is given, 
up to rotations O and O', in the form A = (1 —p)I+p'R. 
It is easily seen that one of the singular values of such 
a transformation is equal to 1 irrespectively of p, while 
the remaining two are identical and given by 1 — 2p. The 




(d) 

FIG. 3. The absolute values of density matrix elements 
(left) and ellipsoids representing the transformed Bloch sphere 
(right) measured by ancilla-assisted quantum process tomog- 
raphy for quantum channels obtained by switching between 
the identity and a selected 180° rotation; (a) the input den- 
sity matrix and the initial Bloch sphere, with selected parallels 
and meridians marked for reference; (b-d) quantum channels 
obtained by changing the duty cycle of the voltage driving 
the controller. Note that the points representing the poles of 
the original sphere do not coincide with the principal axes of 
the ellipsoids, marked with thick straight black lines. 



voltages driving the EPC were programmed in a tem- 
poral loop with a 150 ms period, and the parameter p 
defined the duty cycle during which the EPC operated 
as R. 

The experimental results are shown in Fig. [21[b-d). 
The left column depicts output density matrices recon- 
structed using the maximum likelihood algorithm, while 
the right column represents measured quantum channels 
with the help of ellipsoids that are obtained by acting 
with the channels on the input Bloch sphere. It is worth 
noting that the principal axes of all the ellipsoids are 
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FIG. 4. A graph of the singular values for channels of the form 
(1 — p)I + pR as a function of the duty cycle p. The three sin- 
gular values corresponding to an experimentally characterized 
channel are represented by a square, a circle, and a triangle. 
Solid lines depict values predicted directly from the value of 
the duty cycle as (1,1 — 2p, 1 — 2p) 



oriented along the same direction to within 3°. This is 
consistent with our expectations: the transformation R 
realized by the EPC is applied along the same direction 
determined by the orientation of the mechanical squeezer 
in the EPC, and the birefringence of the input and output 
fibers, which causes rotations of the SOP, is also fixed. 

By means of simple process tomography we found that 
after taking off the voltage implementing R the EPC 
did not return to the state corresponding to the identity 
channel. We attribute this effect to the residual hystere- 
sis in the magnetic elements of the EPC. We partially 
compensated it when implementing channels (1— p)I-FpR 
by setting a small voltage with an opposite sign during 
time intervals when the identity transformation I was to 
be realized. 

The ellipsoids depicted in Fig. [3] contain complete in- 
formation about quantum channels, including unitary 
transformations that accompany depolarization. In or- 
der to analyze more closely the depolarizing characteris- 
tics of the realized channels, in Fig. U] we plot the triplets 
of the singular values versus the programmed duty cycle 
p which defines the share of R in the channel. A good 
agreement between the predicted and measured singular 
values is clearly seen. 

In the second series of measurements, we implemented 
quantum channels that are described by mixtures of the 
identity and two 180° rotations about orthogonal axes, 
which we will denote as Ri and R2 . Such a set of unitary 
transformations allows one to produce all of the trans- 
formations represented by points located on a face of the 
tetrahedron shown in Fig. [1] Using simple process to- 
mography we searched for voltage settings that would 
implement such a pair of rotations. The best approxima- 
tion we were able to find required simultaneous driving of 
all four squeezers. The angle between the rotation axes 
was determined to be 86° ± 3°, with the rotation angles 
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FIG. 5. The absolute values of density matrix elements (left) 
and ellipsoids (right), representing the transformed Bloch 
sphere shown in Fig. [31 a), measured by ancilla-assisted quan- 
tum process tomography for quantum channels obtained by 
switching between the identity and two 180° rotations with 
approximately orthogonal axes. Straight thick black lines 
mark the orientation of the shortest principal axis of an ellip- 
soid. 



given by 175° ± 3°. These findings strongly suggest that 
the alignment of the squeezers inside the EPS must have 
been significantly different from the 45° relative rotations 
that give parametrization of the polarization transforma- 
tion in terms of the Euler angles. 

We used the rotations Ri and R2 to realize a series of 
maps of the form A = (1 - p)I -I- p(Ri -I- R2)/2. Such 
maps are represented by points located on the median 
of one of the faces of the tetrahedron depicted in Fig. [1] 
Analogously to the first series of measurements, the set- 
tings of the EPC were updated in a temporal loop with 
a period of 50 ms. Within a single interval, each of the 
rotations Ri and R2 was realized during a fraction p/2 
of the interval, followed by an application the identity 
I, with voltage settings taking into account the partial 
compensation for the hysteresis. 

Exemplary results are presented in Fig. [5] depicting 
output density matrices and Bloch sphere transforma- 
tions. In order to analyze the results more closely, we 
extracted the singular values and projected them onto 
the tetrahedron face, shown in Fig. [HI Although the ro- 
tations Ri and R2 were applied over identical fractions 
of the cycle, the reconstructed maps do not lie on the 
median of the tetrahedron face. This difference cannot 
be fully explained by the facts that the rotations Ri and 
R2 were characterized by angles deviating from 180° and 
that their axes were not strictly perpendicular. We at- 
tribute the observed disagreement to the hysteresis in the 
squeezing elements of the EPC. Nevertheless, it is worth 
noting that within the experimental errors the triplets 
of the singular values do form a straight line located on 
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FIG. 6. An orthogonal projection of triplets of singular values 
for experimentally implemented channels of the form (1— p)I+ 
p(Ri -|-R2)/2 onto the tetrahedron surface The solid line is the 
least square fit to the experimental points, compared to the 
median marked with a dashed line. The points labelled with 
(a) and (b) correspond to channels depicted in Fig. [5] The 
singular values are given explicitly in the upper left corner of 
the plot. 

a face of the tetrahedron. The mismatch between the 
predicted and the observed characteristics of the channel 
should not pose a substantial difficulty in practice, be- 
cause the averaged operation of the EPC over the entire 



period is repeatable to a very good degree. Therefore 
it is possible to determine, via process tomography, volt- 
age settings that realize a required quantum channel, and 
use them in further experiments. JVIoreover, further im- 
provement could be made by using a larger number of 
squeezing elements, equal to 5 or 6, which, as suggested 
in [l9| . should significantly improve the performance of 
the polarization controller. 

V. CONCLUSIONS 

We used an electrically driven polarization controller 
to implement unital quantum channels acting on the po- 
larization state of a single photon. The characteristics of 
the channels were verified with the help of ancilla-assisted 
quantum process tomography. Although the controller 
requires a calibration procedure for voltage settings to 
implement a given channel, its operation is repeatable to 
a very good degree. Therefore it can be reliably used 
in more complex experiments testings methods to man- 
age decoherence effects in emerging quantum-enhanced 
technologies. Its advantages are low cost, uncomplicated 
driving, and compatibility with fiber-optic setups. 

This work was carried out at the National Labora- 
tory of Atomic, JVIolecular and Optical Physics (FAIVIO) 
in Toruh and was supported by the European Commis- 
sion under the Integrated Project Qubit Applications 
(QAP) funded by the 1ST directorate as Contract Num- 
ber 015848, Polish MNiSZ grant 1 P03B Oil 29, and 
AFOSR under grant number FA8655-06-1-3062. 
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